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Abstract 

In this paper we present a generalized Cassi-Baxter equation to take into account the effect 
of drop pressure on the apparent contact angle 9 app . Also we determine the limiting pressure 
pw which causes the impalement transition to the Wenzel state and the pull-off pressure p ou t 
at which the drop detaches from the substrate. The calculations have been carried out for axial- 
symmetric pillars of three different shapes: conical, hemispherical topped and flat topped cylindrical 
pillars. Calculations show that, assuming the same pillar spacing, conical pillars may be more 
incline to undergo an impalement transition to the Wenzel state, but, on the other hand, they are 
characterized by a vanishing pull-off pressure which causes the drop not to adhere to the substrate 
and therefore to detach very easily. We infer that this property should strongly reduce the contact 



angle hysteresis as experimentally osberved in Ref. [21[. It is possible to combine large resistance to 
impalement transition (i.e. large value ofpw) and small (or even vanishing) detaching pressure p ou t 
by employing cylindrical pillars with conical tips. We also show that depending on the particular 
pillar geometry, the effect of drop pressure on the apparent contact angle 6 app may be more or less 
significant. In particular we show that in case of conical pillars increasing the drop pressure causes 
a significant decrease of 9 app in agreement with some experimental investigations [161 ] . whereas app 
slightly increases for hemispherical or flat topped cylindrical pillars. 

PACS numbers: 47.55.dr, 68.08.Bc, 46.55.+d, 68.08.-p, 65.40.gp, 67.30.hp, 
Keywords: Superhydrorepellence, wettability, lotus effect, microstructured surfaces. 
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I. INTRODUCTION 



Roughness-induced hydrophobicity is a well known effect observed in many plant leaves, 
e.g. Sacred Lotus leaves (Nelumbo nucifera) and in other biological system as water 
striders H] (Gerris remigis) or mosquito (Culex pipiens) eyes \12\- In all such cases the 
surface asperities make the liquid able to be suspended on the asperity tips, resulting in 
a very large contact angle (CA). The superlative water- repellency of such natural surfaces 
would be very appreciated in many micro- and macro- engineering applications, as liquid 
drops on super-hydrophobic surfaces may be very easily moved from one position to the 
other by simply applying an external force, resulting in the possibility to create a chemical 
microreactor and microfluidic microchips [13], 42|. Beside these high-tech applications, there 



is a strong interest of engineers in developing new commercial products as self-cleaning paints 
and glass windows j4j, or super-hydrophobic optically transparent self-cleaning surfaces j9[, 



34j , [351 ] . which may be used as coatings in the automotive field, as car-windshields and 
biker helmets, where impacting rain drops must be easily repelled. Some studies have 
shown, indeed, that falling drops may fully rebound on such water-repellent surfaces with 
very high restitution coefficients ~ 0.9 |40j. As a consequence of the high technological and 
commercial impact of super-hydrophobic surfaces, in the last decade a great deal of research 
has been spent trying to mimic the super-hydrorepellent surfaces of living organisms. Thus, 
many artificial surfaces have been prepared attempting to achieve this objective. Fig. [T] 
shows some examples of super-hydrorepellent man-made surfaces. At a first sight, different 
surfaces may appear equally good candidates to mimic the super-hydrophobic properties 
of Natural surfaces. Indeed, being inspired by the solutions offered by Nature, several 
geometries have been explored, e.g. uniform arrays of flat topped cylindrical pillars 0] [Fig. 
Ufa)], or tapered asperities [2l(, [see Fig. GJb)], or even nanotube forests with rounded 
tips 20] as shown in Fig. [He). Beside these few examples, also super-hydrophobic fractal 
surfaces have been produced [lOj- 41] . which show apparent contact angles (CAs) up to 174°. 
However, the different shapes of such microstructured surface, which reflect the variety of 
Natural solutions, should also have some practical implications which makes them not really 
completely equivalent. Liquid drops on super-hydrophobic microstructured surfaces ma y b e 



observed mainly in two different states (although intermediate states may also exists [22| 
[l|) the Cassie-Baxter 0] and the Wenzel j43| states. A drop in a Cassie-Baxter state is just 
suspended on the asperities of the underlying surface, which therefore behaves as a fakir- 
carpet. The Wenzel state is, instead, characterized by complete contact between the drop 
and the substrate. The Wenzel state is usually unwanted as it results in a strong adherence 
between the drop and the substrate and in a very pronounced CA hysteresis 16], 3{|. Thus, 
in order to prevent strong adhesion between the drop and substrate one has to design the 
super-hydrophobic microstructured surface in such a way to prevent the Wenzel state to 
be formed and make stable the fakir-droplet state. In such sense a first attempt was made 



in 



231 ] where a model employing the capillary rise of a liquid in contact with a stripwise 



heterogeneous surface was developed to study the effect on contact angles. More recently 
in [3] two criteria were proposed to compare the superhydrorepellent properties of different 
microstuctures from a wetting point of view. However, the proposed criteria do not take 
into account the effect of the internal pressure of the droplet, which is strictly related to 
surface tension and curvature of the air-liquid spherical interface and therefore to its volume. 
An interesting new approach with molecular dynamics was proposed in to study the 
behavior of liquid nanodropltes on rough surfaces. However, as shown by different authors 
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FIG. 1: A few examples of artifical superhydrorepellent surfaces, (a) A uniform array of very 
slender cylindrical pillars with a flat tip (adapted from Ref. [5]), (b) a uniform array of cylinders 



carbon-nanotube forest (adapted from Ref. 20]) 



with tapered shaped tips (adapted form Ref. [2l|]), and (c) an example of superhydrorepellent 



16 . S], (26l |. [39 ] . the drop pressure may have a critical role in determining if composite 
interface may be formed at the interface between the drop and the microstructured substrate. 
Large drop pressures may be generated during the impact of drops on the substrate, and 
in this shown in Ref. [40j, |25J, [26( high impact velocities (i.e. large impact 



pressures) may destabilize the fakir-state, cause the transition to the Wenzel state and 
make the droplet not able to bounce. Fig. [2] shows indeed two water drops of the same 
volume on a microstructured super-hydrorepellent surface. The drop on left has been gently 
deposited on the substrate, thus allowing the formation of a fakir-droplet state, the drop 
on the right has instead undergone a transition to the Wenzel state as a consequence of an 
increase of the drop pressure above a critical threshold. 

The transition between the composite (Cassie) and wetted (Wenzel) states has been 



3, @, H, H, H), but usually the 



investigated theoretically in many papers ([271]. 
effect of the liquid drop pressure is not taken into account. Only very recently an interesting 
study [33[ has been presented where the superhydrorepellent properties of a surface with 
cavities has been investigated and the effective energy of such systems has been studied by 
also including the influence of drop pressure. In this paper we study the system by means of 
an energy approach of which the most general formulation has been given by Lipowsky [l7l ]. 
who by means of a minimization technique proposed a generalization of the Cassie-Baxter 
and Wenzel laws for liquid drop sitting on chemically heterogeneous but flat substrates, 
where the area fraction of each phase was given a priori. Also, there are others paper treat ing 



45 



the pro blem of hydraulic pressure in determining the stability of Cassie-Baxter state 
and [l9|], for example, treated this problem in the case of pillars with flat tip with sharp edges 
so that the fraction of the projected area that is wet is assigned a priori. In particular, [45] ] 
gave an expression of the critical hydraulic pressure for which the transition process between 
the Cassie-Baxter wetting mode and the Wenzel one occurs, 19| also studied the effect of the 
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FIG. 2: Two millimetric water drops of the same volume on a microstructured super-hydrorepellent 
surface (adapted from Ref. [a]). A fakir state (i.e. Cassie-Baxter state) is obtained by gently 
placing the drop on the substrate (on the left). A full-contact condition (Wenzel state) is obtained 
by incresing the drop pressure, e.g. by let the drop impacing the surface with enough speed (on 
the right). 



contact line length by varying the shape of the pillars while maintaining constant the area. 
Indeed for rounded pillars the wetted area is not fixed a priori and should be determined 



as a part of the problem through a minimization technique. In [18] the authors studied the 
case of a surface with an array of small spheres on it and determined the pressure required 
to force the meniscus to bend sufficiently to touch the underlying surface. However, if the 
pillars are sufficiently tall the transition to Wenzel state may be achieved before the meniscus 
touches the underlying surface because of thermodynamic instability. 

In Ref. j§| it has been shown that to prevent the impalement transition it is necessary 
to increase the critical pressure pw at which this transition occurs (we refer to pw as the 
Wenzel pressure); high values of pw are indeed a strict requirement in those engineering 
applications where rain falling drops have to be supported by the substrate. However, in 
many cases we also want the liquid drops be very easily detached from the substrate, i.e. 
the pull-off pressure p out be reduced almost to zero. This property is, for example, found 
in many insects, that usually walk or skate on the free liquid-air surface. Such insects not 
only need not to sink into the water, but also need to detach easily from the liquid free and 
surface move and run easily on it. 

In a preceding paper jl] one of the authors has analyzed the wetting/non- wetting prop- 
erties of a liquid drop in contact with an extremely idealized ID rough profile, i.e. a simple 
sinusoidal profile. The analysis have clarified some theoretical points (mainly from a qual- 
itative point of view) of wetting non wetting behavior of super-hydrorepellent surfaces. 
Here the study is extended to 2D microstructured surfaces with periodic a distribution of 
axial-symmetric micropillars, which on the other hand are very commonly utilized in such 
applications. Some hints to design such super-hydrophobic surfaces to achieve both the aim 
of large pw values and low pull-off pressures p out are also provided. 
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FIG. 3: Geometry of the conical (a), emispherical topped (b), and cylindrical (c), pillars. 



II. FORMULATION 



Let us consider a periodic distribution of chemically hydrophobic (e.g. fluorinated) pillars 
(thermodynamic contact angle 9 e > 7r/2). We assume that the elementary cell of the periodic 
structure is a square, although we can deal with any type of periodic distribution of micro- 
pillars. We develop the analysis for three types of micro-pillars: conical pillars, hemispherical 
topped cylindrical pillars, and flat topped cylindrical pillars. We assume that the micropillars 
are stiff enough to consider negligible their deformation under the action of the drop pressure. 
This, indeed, is a very good approximation in many practical cases, and is always employed 
in theoretical investigations dealing with super-hydrorepellence. We also assume that the 
liquid is incompressible (i.e. we neglect the contribution of the liquid elastic energy) and 
the drop is very slowly evaporating (i.e. the time scale to reach the equilibrium is much 
shorter than the time scale of the evaporation). We observe that, being the diameter of the 
liquid drop in the range of millimeters, whereas the linear spacing 2A between the pillars 
in the range of micro- or even nano-meters, the drop can be considered as a semi-infinite 
liquid space when we analyze the problem at the microscale. Thus, when we look the 
drop-substrate interface at very large magnifications the state of the systems is completely 
determined by the following state parameters: the drop-pressure p at the liquid-substrate 
interface, the real solid-liquid contact area, the liquid-air free surface at the interface, and 
the penetration A (see below) of the liquid drop inside the pillar forest. 

Fig. [3] shows the geometry of the pillars and the parameters we use to describe the 
position of the triple line. The reference plane (x, y, 0) is placed at the base of the pillars, 
and the z-coordinate is directed toward the top of the pillars. The ^-coordinate of the free 
liquid surface will be referred to as u (x, y). Because of the substrate corrugation, the liquid 
can either wet the whole substrate surface or be in stable or metastable partial contact 
with it. In case of partial contact the liquid/air interface must satisfy the Laplace formula. 
Assuming that the slope of the liquid-air interface is sufficiently small, which simply requires 
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that p\j {2^ la) "C 1, the Laplace formula reads 



V 2 w (x, y) = u xx (x, y) + u yy (x, y) 



P 



(1) 



lLA 



where u xx = d 2 u/dx 2 , u yy = d 2 u/dy 2 and jla is the liquid-air surface tension. The condition 
p\/ (Z^la) <C 1 is satisfied in most cases, e.g. in case of water 7^ = 72 mJ/m 2 and assuming 
A ~ 1/ini one obtains p < 1.4 bar which is in most cases true. This then implies the angle 
which the free surface forms with the pillar, for the geometry we have investigated and 
R < A (where R is the radius of the pillars), need to be less than 30°. 

Eq. ([1]), because of periodicity, can be solved over a quarter of the elementary square cell. 
However we need also to specify boundary conditions. A first boundary condition has to be 
written at the triple-line, which represents the contour delimiting the liquid-solid interface. 
Let us call the projection of this contour on the (x, y, 0) reference plane with the symbol 
L. We observe that, in general, the curve L (for which the mathematical expression in 
implicit form can be given as fi (x, y) = 0) is not known a priori and has to be determined 
by requiring that the total energy of the system is stationary at equilibrium (see below). 
Therefore at the triple line the following equation must hold true 



where L = {(x, y) G 9? 2 | /l (x, y) = 0}, and h (x, y) is the function describing the shape of 
the pillars. Eq. ([2]) simply states that the liquid-air interface and the solid-liquid interface 
must intersect at the triple line. Also the following Neumann boundary conditions must be 
satisfied to account for periodic conditions 



where u x = du/dx, u y = du/dy, x satisfy the condition /l(x ,0) = and, similarly, y 
satisfies the condition f L (0,y ) = 0. 

Observe, that for flat topped cylindrical pillars the function (x, y) is known a priori 
being simply /x (x, y) = x 2 + y 2 — R 2 , where R is the radius of the pillar. In the other two 
cases fi {x, y) is not known a priori and must be determined as a part of the solution of 
the problem. Indeed, the physical problem, we are dealing with, belongs to the class of free 
boundary problems, and requires an additional condition to achieve the complete solution. 
This additional condition is simply the requirement that at equilibrium, for any given drop 
pressure p, the total energy of the system (in our case the Gibbs energy G) is stationary. 
Of course, in the general case, the Gibbs energy is a functional defined on the vector space 
of functions fi(x,y), and one should require that it is stationary at equilibrium to find 
the Euler-Lagrange equations and determine the quantity y), making the problem 

belonging to the class of variational problems. Therefore, the complete solution seems to be 
very complicated and expensive from a numerical point of view. However we can strongly 
reduce the complexity of the problem if we recall that p\/ (Z'Jla) <C 1, in such a case the 
slope of the free liquid-air surface is small. This implies that also the slope of the contour 
representing the triple line is small and we conclude that under this assumption the triple 
line will only negligibly deviate from a circumference in case of axial-symmetric pillars ( [19| | 
showed the variation of the angle that the free surface forms with the pillars, due to the 
tortuosity of the triple line, is negligible (being limited to only 3%). Thus, the unknown 



u{x,y) = h (x,y) ; 



(x,y) e L 



(2) 



u x (0, y>y ) = 0; u x (A, y) = 
u y (x > xq, 0) = 0; u y (x, A) = 



(3) 
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FIG. 4: The domain of integration of Eq. (fT]). 



function (x, y) takes a much simpler form f L (x, y) = x 2 + y 2 — r%, where r is the unknown 
radial position of the triple-line, and the total energy of the system simply becomes a function 
of the free parameter r^. Thus, ro can be determined by enforcing the equilibrium condition, 
i.e. requiring that dG/dr Q = 0. However, despite the apparent simplicity of Eq. (pQ), the 
particular shape of the domain of integration (see Fig. HJ) does not make it possible to obtain 
a solution in closed form. We, therefore, have employed a finite difference approach to solve 
Eq. (|T]) with the mixed boundary conditions Eqs. (j2J) and (|2D. 

For isothermal conditions and constant pressure the Gibbs energy is 



where A (r ,p) is the penetration of the rigid substrate into the semi-infinite liquid and H 
is the Helmholtz free energy. The penetration A is given by 



where s is the average profile of the liquid. Eq. (j^J) represents the equation of state of the 
system since it allows to determine one of the three quantities A, ro, p once known the other 
two. To define the Gibbs energy we first need to express the Helmholtz free energy H which in 
our case is just the total surface energy of the system H (r ,p) = •JlsSls + IlaSla + IsaSsa, 
where 'JlsSls: IlaSla and JsaSsa are the surface energies at the interfaces liquid/solid, 
liquid/air and solid/air, respectively. Utilizing the Young's equation jLA C0S ^e J rlLS~ ISA — 
0, with 9 e being the Young's CA at equilibrium, the Helmholtz free energy reads 



G (r , p) = H (r , p) - p\ 2 A (r , p) 



(4) 



A (r ,p) = h-s (r ,p) 



(5) 



H (r , p) = Jla (Sla ~ S LS cos 9 e ) + j S aS s 



(6) 



and the Gibbs energy becomes 



G (r , p) = Jla {Sla - S LS cos 9 e ) - p\ 



2 A + jsaSs 



(7) 



In Eqs. (EJ) and (|7j) Ss represent the total surface of the substrate over the single square 
cell. We now require that at fixed load (i.e. at fixed drop pressure p) the Gibbs energy 
is stationary to enforce equilibrium conditions and close the system of equations with the 
following condition 

£ = (8, 

Stability or instability of equilibrium can be easily determined by looking at the sign of 
d 2 Gjdr 1 : Local stability is guaranteed when the energy has a local minimum, i.e. when 
d 2 G/dr 2 > 0, whereas instability is detected when d 2 G/dr 2 < 0. 



A. The apparent contact angle 

At the macroscopic scale the miscrostructure of the substrate is not observable, therefore 
at the macro-scale the observer will measure an apparent liquid solid surface energy { r )Ls) e ff 
equal to 

{lLs) eff = H/X 2 (9) 

using Eq. © (lLs) eff becomes 

/ \ 1la(S L a- S ls cos 9) 

VYLS)eff = p + USA) eff (10) 

where we have defined (jsa)^* = (jsaSs) as the effective solid-air interfacial energy. 
At equilibrium the above definition Eq. (1101) allows to write a modified Young's equation 
and evaluate the apparent contact angle 9 app as 

_ (7ls) £// - (lSA) eff _ S LA -S LS cos6 

COS U a pp — — l-L-U 

ila y 

The above Eq. ffTTj) represents a generalization of the Cassie-Baxter equation [7], which 
takes into account the influence of the interfacial drop pressure, and, of course, holds true 
only at equilibrium. 

The above simple considerations make clear that a design criteria based on the maximiza- 
tion of the apparent contact angle is effective only if we include the effect of the pressure. 
Therefore, an optimization of the surface topography only based on the apparent contact 
angle evaluated by the Cassie-Baxter equation can lead to misleading conclusions. 



III. EQUILIBRIUM CONDITION: A SIMPLIFIED APPROACH 

The general approach considered above which leads to an optimization procedure to find 
Gibbs energy stationary points is very time consuming. For this reason here we present 
a much simpler procedure to determine the equilibrium conditions without the need of a 
complicated and time-consuming minimization procedure. Indeed, if we were able to enforce 
the equilibrium condition before solving Eq. (pQ), we would have the advantage of strongly 
reducing the complexity of the problem, since, as we are going to show, this allows to 
readily calculate the contact radius r for any given drop pressure p and hence to 'bypass' 
the minimization of the Gibbs energy. This, indeed, can be done by requiring that the CA 
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FIG. 5: The forces acting on the emispherical topped pillars, (a); and on the conical pillar, (b). 
Note that we have assumed that the triple line contact is circumference of radius r$ and that the 
CA angle is just equal to the Young's contact angle 9 e . 

at the triple line is just equal to the Young's contact angle 9 e . Therefore, as shown above, 
that the triple line negligibly deviates from a circumference this condition allows to calculate 
the real liquid-solid contact area for any given drop pressure. Of course the determination 
of the complete thermodynamic state of the system still requires the calculation of the 
free liquid-air interface u (x, y) to determine the total interfacial energy and therefore the 
apparent contact angle 9 e ff and penetration A. However, this calculation can be carried 
out a posteriori by solving Eq. (jTj), with the conditions ([2]) and ([3]), without handling the 
minimization problem of the total energy. 

A. Conical pillars 

Let us observe Fig. (b) where the forces acting on a conical pillar are shown to be the 
force F that the rigid substrate applies to the pillar, the surface tension ^la at the triple 
line, and the pressure p of the liquid. Thus, the equilibrium writes as 



where a is the half-cone angle. Now let be A the measure of the area covered by the 
elementary cell (not necessarily square) of our periodic distribution of asperities, (in the 
case of a square cell the quantity A is simply A = A 1 / 2 / 2), and enforce the equilibrium of 
the whole single cell. Because of the Neumann boundary condition Eq. (J^D the liquid-air 
surface tension at the outer boundaries of the single cell will not give any contribution to 
the equilibrium along the z-direction. Therefore we can write F = pA and the above Eq. 




(12) 



(IT2"|) becomes 



7r (r/2) cos (6*6 — a) 



(13) 



p = 



1 -7r(f/2) 2 



where f = ro/A and p = pX/'jLA, being A = A 1 / 2 /2 a characteristic length. Notice when 
a > 9 e — 7r/2 the pressure becomes negative and the drop wets the substrate in a Wenzel 
state. Hence the angle a has to satisfy the condition < a < 6 e — 7r/2 to guarantee a 
positive value of the drop pressure p > 0. In such case the drop pressure p continuously 
increases with the radius f , i.e. stability of equilibrium is always guaranteed. The pull-off 
pressure can be easily evaluated as p out — p (r — 0), which as shown by Eq. ([TBI is zero, i.e. 
drops on conical pillars can be very easily detached from the substrate. We observe that 
in all practical cases the conical tip will never present a real sharp corner. However this 
does not change our conclusion. The conical pillar will always guarantee a negligible pull-off 
pressure, as the radius of curvature of its rounded tip is always much smaller than radius of 
the pillar itself. Moreover, notice that since the effect of the real extension of the liquid-solid 
contact area on adhesion between the drop and surface is already taken into account when 
we calculate the contribution of the surface energy to the total energy of the system the 
detachment of a liquid drop from a conical pillar requires a vanishing small force even if the 
liquid-solid contact area can be larger than in cases of hemispherical or flat topped pillars. 

B. Hemispherical topped cylindrical pillars 

For hemispherical pillars (see Fig. [5] (a)), we can follow the same procedure as outlined 
above to write the equilibrium 



where R is the radius of the sphere and ip represents the angular coordinate of the liquid- 
pillar triple line [see Fig. [51(a)] . Eq. (1141) can be conveniently rewritten in a dimensionless 
form as 



where we have still used that p = pX/^LA and R = R/X. At fixed load, stability implies 
dp/ dip > which is equivalent to have d 2 G/dr 2 l > 0, whereas unstable equilibrium conditions 
are characterized by value of dp/ dip < 0, i.e. d 2 G/dr 2 l < 0. Therefore, the threshold value of 
pressure at which instability occurs can be determined by enforcing the condition dp/ dip = 0, 
i.e. 

— ttR 2 sin 9 e — 8 cos (2<p>) sin 9 e + ttR 2 cos (2(p) sin 6*6 — 8 cos 9 e sin (2ip) = (16) 

We observe that two unstable conditions can be found in general: the first one corresponds 
to the impalement transition to the Wenzel state and is reached when the pressure increase 
over a value pw (see below), the second one is achieved when the pressure decreases below 
a value p out which we call pull-off pressure. In particular we stress that the pull-off pressure 
p ou t represents the lowest value of the pressure at which it is still possible to find a stable 
minimum of the total energy of the system. If the liquid pressure decreases below this 
threshold value the liquid drop detaches from the substrate. However the drop does not 
detach as a whole, but rather via the edge propagation of the triple line toward the inner 



27r (i?sin ip) sin (8 e + <p) + F — n (i?sin ip) 2 p = 



(14) 



V 




(15) 



1 - (tt/4) ^sinV 
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FIG. 6: The two limiting states that a liquid drop in contact with a flat-topped cylindrical pillars 
can assume. Limiting conditions at pull-off, (a); and limiting condition at the transition from 
Cassie-Baxter to Wenzel state, (b). 



region of the liquid-pillar. We also observe that from a conceptual point of view, the pull-off 
pressure defined in the present paper is analogous to the maximum detachment force defined 



in Ref. 24] for capillary bridges. In fact, in our case the thermodynamic contact angle is 



9 e > tc/2 and in this case the maximum detaching force for capillary bridges occurs exactly 



at the transition I point defined in Ref. 24| (see figure 5 of the cited paper). Beyond this 
point stable contact conditions cannot be achieved and the drop will necessarily detach from 
the substrate. 

As in case of hemispherical topped pillars < p < ir/2, the equation dp/ dip = gives 
just one solution (p = ip out , which corresponds to the pull-off condition p out = p ((p out )- The 
impalement transition to the Wenzel state is instead obtained at (p — (fw = it/ 2 

Pw =PVPw) = -cos9 e - T — (17) 

1 - Ti R 2 A 



C. Flat topped cylindrical pillars 



In case of flat topped cylindrical pillars the liquid-solid contact area on each pillar is 
fixed, and equal to ttR 2 . In such case, depending on the drop pressure values, the slope of 
the liquid profile at the triple line can vary between two different limiting values as clearly 
shown in Fig. EJ We can, therefore, easy determine the critical pressure p out at pull-off [see 
Fig. EUa)] as 

„ kR/2 . . 

Pout = - sin 6 e 7—- (18) 

1 - txR 2 /A 
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and the critical dimensionless pressure pw [see Fig. MJo)] as (see also [22 ]) 

n 7T-R/2 . . 

p w = - cos e - ^— (19) 

1 - nR?/A 



which, as expected, is exactly equal to the value found for hemispherical pillars [see Eq. ([Tj 
It is noteworthy that, in case of flat topped cylindrical pillars, the ratio \p 0U t/Pw\ — tan^g 
non depending on the cylinder radius. 



IV. RESULTS 

In what follows we assume that the Young's contact angle is 9 e = 109°, i.e. we assume 
that the underlying microstructured surface is chemically hydrophobic. 



A. Conical pillars 

Conical pillars have been analyzed assuming that R = R/X = 0.5, where R is the radius 
of the base circle. Fig. [7J shows the contact radius f as a function of the dimensionless 
drop pressure p for different values of the dimensionless pillar height h. The black solid 
curves are obtained by Eq. ffTBl) whereas the red-dashed ones are obtained by minimizing 
the total energy of the system as explained in Sec. [Til Notice the very good agreement 
between the two approaches. Also Fig. [JJ shows that f initially increase proportionally 
to the drop-pressure p, as indeed predicted by Eq. (TT3|) . however, as the drop pressure is 
increased further, the contact radius also increases and the denominator in Eq. (TT3|) may 
become not negligibly smaller than one. This, in turn, causes the drop pressure p to increase 
more than linearly with f thus explaining the deviation of curves in Fig. [7J from linearity. 
Since for conical pillars the partial fakir-droplet state is always stable (i.e. increasing the 
pressure does not force the system to undergo a spontaneous transition to the Wenzel state), 
a different threshold pressure pi has to be defined. This is simply the drop pressure at which 
the free liquid-air interface touches for the first time the basal of the pillars. This condition, 
indeed, has been shown experimentally to easily trigger a sharp transition to the Wenzel 
state [22|. Fig. [8] shows the dimensionless limiting pressure pL as a function of the pillar 
aspect ratio h. Observe that increasing h also increases the limiting pressure pi, since (at 
fixed A) the liquid-air interface will be farther from the bottom of the pillar. However, p^ 
cannot increase above the limiting asymptotic values obtained for h — > oo, i.e. a — > 0. 
When this happens, the conical pillar becomes an infinitely tall cylinder of dimensionless 
radius R for which the Cassie-Baxter state becomes unstable when the drop-pressure reaches 
the values pw given by Eq. (TT9|) which indeed is just the asymptotic value (pz,)^ shown by 
the dashed line in Fig. [HJ 

Figure M shows the dimensionless penetration A = A/ A of the liquid drop into the pillar 
forest as a function of the dimensionless pressure p. As expected A increases as the pressure 
p and aspect ratio h are increased. Indeed, increasing h, at fixed R and A, makes the 
cone sharper and sharper. This would in turn reduces the liquid-solid contact radius if the 
penetration were maintained fixed. But, since the dimensionless pressure p at equilibrium 
decreases if the contact radius f is reduced [see Eq. (TT5|) ]. the drop actually needs to increase 
the penetration to increase f and, thus, sustain the applied pressure p. 



12 









0.1 0.2 
P 



0.3 



FIG. 7: The dimensionless radius of contact r as a function of the dimensionless pressure p at 
equilibrium, for a drop in contact with a periodic distribution of conical pillars. The red-dashed 
curves are results of the minimization procedure as explainen in Sec. II, black continuous lines 
are instead determined by the simplified approach of Sec. III. The matching between the two 
methodologies is excellent. Result are presented for different values of h = 2, 4, 6, 8, and R = 0.5, 
the Young's contact angle is 6 e = 109°. 



It is noteworthy to observe that if the dimensionless height of the conical pillar is strongly 
reduced, we should expect an increase, instead of a decrease, of the dimensionless penetration 
A as a function of h. Indeed, in such case the half-cone angle a would increase toward the 
limiting value 6 e — ir/2 at which the drop spontaneously undergoes a transition to the Wenzel 
state (which is obviously characterized by A = h). This increase of the penetration, as a 
consequence of strong reduction of h, is, indeed, also observed in Fig. [9] for h = 2, where 
the corresponding A values are larger than those obtained for h = 4 over a large range of p. 

Figure [TD] shows the calculated apparent contact angle 9 app [see Eq. (TTTT) ] . as a function of 
the radius of the base circle R and different drop pressures p. In the limit case of zero pressure 
apP = 180° independently of R, since in this case the stable state is a perfect Cassi-Baxter 
state with the drop just touching the tip of the conical pillars. Increasing p the apparent 



contact angle decreases in qualitative agreement with some experimental observations [16 
The reason of this decrement is related to the significant increase of penetration A which 
in turn determines a strongly increase of the energy term pA a therefore a reduction of the 
apparent contact angle [see Eq. fllip ]. Notice that for each given pressure p > a minimum 
value of the base radius is needed in order to stabilize the composite interface and avoid the 
transition to the Wenzel state. This explains why the curves at different drop pressures in 



13 



0.3 



0.2 



0.1 











(M 



R =0.5 

0e = 109° 



10 



20 



30 



h 



FIG. 8: The threshold pressure pl as a function of the pillar aspect ratio h. Results are presented 
for R = R/X = 0.5 and 9 e = 109°. Observe that increasing h the limiting pressure reaches an 
asymptotical value {pl)^ which is just the values given by Eq. ([Tj 



Fig. [10] are plotted starting from different values of the base radius R. 



B. Hemispherical topped pillars 

To analyze the behavior of hemispherical topped pillars we use the simplified approach 
described in Sec. Ill to calculate the area of liquid-solid contact as a function of pressure, 
then we solve Eq. ([T|) to determine the shape of the free liquid-air interface and therefore 
the penetration of the liquid drop and the apparent contact angle. Fig. [TT] shows the di- 
mensionless liquid-pillar contact radius f at equilibrium [Fig. [TT] (a)] and the dimensionless 
penetration A [Fig. [TT1(b)]. as a function of the dimensionless drop pressure p. Two types of 
lines are shown. Solid lines are stable equilibrium branches, whereas dashed lines represent 
unstable branches at fixed load. As expected, also in this case the contact radius and the 
penetration increase with the applied drop pressure, but this time the f vs. ft law strongly 
differs from being linear. Finite negative values of \p ou t\ are related to the finite size of the 
liquid-solid contact that still exists when the drop is about to detach from the substrate. 
This is often strongly unwanted since, beside the large detaching force, it usually leads also 
to large contact angle hysteresis 



15 , 38 



Also observe that for hemispherical topped cylin- 
drical pillars impalement transition occurs spontaneously when the drop pressure reaches 
the limiting value pw (at which the penetration A jumps to h), in contrast with what we 
have found for the conical pillars. Therefore, for hemispherical pillars the height should be 
chosen by taking care that at p — pw the free liquid-air interface does not touch the bottom 
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FIG. 9: The dimensionless drop penetration A, as a function of the drop pressure p. Results are 
presented for R = 0.5 and 6 e = 109°, and for different values of the conical pillars aspect ratio h 
= 2, 4, 6, 8. 
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FIG. 10: The apparent contact angle 6 app as a function of the base radius R for conical pillars. 
Results are shown for different drop pressures p. The aspect ratio of conical pillars is h = 8. 
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FIG. 11: The dimensionless radius r at equilibrium, (a); and the dimensionless penetration A 
at equilibrium, (b); as a function of the dimensionless drop pressure p for hemispherical pillars. 
Full lines are stable branches, wherease dashed lines are unstable branches at fixed load (at fixed 
penetration the stability extends till to the point where A has a minimum). Results are presented 
for different values of the dimensionless sphere radius R = 0.2, 0.3, 0.4, 0.5 and for 8 e = 109°. Each 
curve ends at a certain value of p = pw, which only depends on the pillar radius and contact angle 
9 e . When this value of drop pressure is reached the penetration sharply jumps to the unit value, 
i.e. the drop ungergoes a sharp transition to the Wenzel state. 

of the pillars forest. 

In case of hemispherical topped pillars, Fig. shows that the apparent contact angle 
slightly increases with pressure p because in this case the liquid/air surface term JlaSla 
increases more than the other two energy terms p\ 2 A and JlsSls- Fig. 02] also shows the 
original Cassie-Baxter solution which corresponds to p — 0. 

C. Flat topped cylindrical pillars 

In this case, for any value of the drop-pressure between the two limits p out given by Eq. 
f|T8|) and pw given by Eq. fTT9|) . the contact liquid-pillar area will be always equal to nR 2 
with R the radius of the cylinder. This makes the problem ([1]) ([2]) and ([3]) linear, and, in 
turn, leads to a direct proportionality between the drop pressure p and the penetration A, as 
indeed confirmed experimentally in Ref. (22|. Notice that, since \p ou t/pw\ = |tan | and 6 e 
is usually not larger than 120°, the pull-off pressure is always larger than p w . Therefore we 
expect that, although a forest of flat topped cylinders can stabilize the Cassi-Baxter state, 
a drop suspended on such a microstructured surface is relatively difficult to detach from it 
and should suffer of strong contact angle hysteresis. 

Fig. (|T3l) shows the variation of the apparent contact angle 6 app as a function of the 
cylinder radius R for different drop pressures. Similarly to the case of hemispherical pillars 
the apparent contact angle grows with the pressure, although in this case this increment is 
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FIG. 12: The apparent contact angle 9 app as a function of the radius R for hemispherical topped 
pillars. Results are shown for different dimensionless drop pressure p. 
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FIG. 13: The apparent contact angle 6 app as a function of the radius R. Results are shown for 
cylindrical pillars and different dimensionless drop pressure p. 



less pronounced. 



V. DISCUSSION AND DESIGN SUGGESTIONS 



Very robust super-hydrorepellent surfaces should possess the ability to support large drop 
pressures, and should also allow the drops to easy abandon the substrate, roll on it with 
almost zero contact angle hysteresis, or even easily bounce on it. As already stated, these 
properties are very desirable in applications such as micro-fluidic-chips and micro-chemical 
reactors, where drops have to be easily moved and positioned. But, also at the macro- 
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scale, they represent a strict requirement for self-cleaning windows, water-hydrorepellent 
windshields, or even hydrorepellent clothing. In these latter cases, the droplet pressure may 
reach large values as a consequence of large inertia forces due to the impact of rain drops. 
As an example, assuming the rain drop falls at the speed of about vq = 10 m/s we can 
easily estimate the maximum impact pressure p max ~ P^o = 1 x 10 5 Pa. Therefore, a self- 
cleaning super-hydrorepellent window should be necessarily characterized by large values of 
the critical Wenzel pressure such that pw > Pm&x- Thus, assuming that the substrate is 
constituted of hemispherical topped cylindrical pillars distribution, and recalling Eq. (fl9j) 
we found 

A < Pw = - cos 6 e 7— = A max (20) 

Pmax 1 - /4 Pmax 

So taking for R the value 0.5, recalling that the liquid-water surface tensions is jla ~ 
72 x 10" 3 J/m 2 , and assuming a thermodynamic contact angle 6 e = 109°, one ends up with 
the value A < 230 nm, i.e. nanotube forests [20j] should be employed in such applications, 
provided that the pillar height is sufficiently large to avoid direct contact between the free 
liquid-air interface and the bottom of the solid surface (we also observe incidentally that 
nanometer spacing between the pillar is also necessary not to alter the transparency of the 
glass). However, we underline that the present paper is focused on steady-state conditions 
and a correct treatment for impactig drops should take in account that the curvature of the 
free surface is determined by the non uniform pressure distribution at the interface. 

In case of conical pillars the critical pressure pi is determined by the condition that the 
air-liquid interface touches the pillar forest ground. It has been observed that the pressure 
pi is always smaller than the critical Wenzel pressure pw calculated for hemispherical and 
flat topped cylindrical pillars, and approaches this value only for infinitely large values of 
h. Therefore, one may be tempted to conclude that if large drop pressures have to be 
supported, the conical pillar shape is not a viable solution. However, we can easily correct 
for this by simply slightly modifying the pillar design to turn it in a cylinder with a conical 
tip on the top, i.e. in a conical topped cylindrical pillar. The new conical pillar design 
guarantees the same critical Wenzel pressure pw as given by Eqs. fTT7|) or ffT9l . but has 
the fundamental benefit to present a vanishing or negligibly pull-off pressure p ou t- As a 
consequence, drop on a forest of conical topped cylindrical pillars should not suffer from 
CA hysteresis, as indeed experimentally observed in Ref. |2l|, and should be able to easily 
roll or slide on the substrate or even to bounce on it with very high restitution coefficients. 
This, should also explain why some biological systems as water striders, which usually walk 
on the free surface of the water, possess a conically shaped distribution of asperities on their 
super-hydrorepellent legs pj| . Also note that it is useless to have pillars taller than the 
minimum height necessary to prevent the contact of the liquid surface with the bottom of 
the pillars forest, in contrast to what is often asserted, i.e. that the taller the pillars the 
more the super-hydrorepellence of the surface. The argument, which is usually provided, is 
that as the pillar height is increased, more energy has to be spent to push the drop in full 
contact with the substrate. Thus, making pillars taller and taller should lead to a very large 
resistance against the impalement transition [38}. However, one should observe that this 
energy can always be provided by the pressure pw acting inside the drop independently of 
the pillar height, and that the real critical condition for the transition to the Wenzel state to 
occur is p = pw- Only in case of very small drops (diameter comparable with the spacing A 
between the pillars) the physical scenario may slightly change. In this case, it may be shown 
that the pillar height can actually play an additional role: because of mass conservation, an 
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additional resistance against impalement transition and drop penetration is generated 
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VI. CONCLUSIONS 

In this paper the behaviour of a liquid drop on super-hydrorepellent surfaces constituted of 
periodic distribution of pillars has been analyzed. In particular the critical drop pressure pw 
which destabilize the fakir-droplet state causing a transition to a Wenzel (full contact) state 
and the critical pressure p out which causes the detachment of the drop from the substrate have 
been studied for three types of periodic micro-structured surfaces: conical, hemispherical 
topped and flat topped cylindrical pillars, regularly disposed on a rigid substrate. Both pw 
and p ou t are equally important to assess the superhydrorepellent properties of surfaces. In 
fact high pw values are requested in all those applications in which very high pressures must 
be supported, e.g. self-cleaning glasses and super-hydrorepellent windshields, whereas small 
values of p out are desirable to guarantee very small contact angle hysteresis and allow the 
drop to easily move on the substrate (e.g. microfluidic chemical reactor, micro-fluidic-chips) 
or, in case of impacting drops, easily rebound from it (e.g. self-cleaning windows, super- 
water- repellent windshields and biker helmet visors). We have shown that the conical pillars 
have a pull-off pressure p out vanishing small (that is an advantage in those applications in 
which liquid drops have to be easily removed from the surface), but the Cassi-Baxter state 
is destabilized for pressures smaller than the critical value pw found for hemispherical or 
flat topped cylindrical pillars. However, a surface microstructured with cylindrical pillars 
with conical tips would have both advantages of large pressure pw and zero pull-off pressure 
p out . Finally, the effect of the pressure on the apparent contact angle 9 app has been studied. 
The analysis has shown that 9 app reduces significantly with the pressure in case of conical 
pillars (in agreement with previous experimental observation), whereas it slightly increases 
for hemispherical or flat topped cylindrical pillars. 
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